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ABSTRACT 

In order to fully understand the gravitational collapse of molecular clouds, the star formation process and 
the evolution of circumstellar disks, these phenomena must be studied in different Galactic environments with 
a range of stellar contents and positions in the Galaxy. The young massive association Cygnus OB2, in the 
Cygnus-X region, is an unique target to study how star formation and the evolution of circumstellar disks 
proceed in the presence of a large number of massive stars. We present a catalog obtained with recent optical 
observations in r, /, z filters with OSIRIS, mounted on the 10.4 m GTC telescope, which is the deepest optical 
catalog of Cyg OB2 to date. 

The catalog consist of 64157 sources down to M = 0.15 Mq at the adopted distance and age of Cyg OB2. 
A total of 38300 sources have good photometry in all three bands. We combined the optical catalog with 
existing X-ray data of this region, in order to define the cluster locus in the optical diagrams. The cluster locus 
in the r - i vs. / - z diagram is compatible with an extinction of the optically selected cluster members in the 
2.64"' < Ay < 5.57'" range. We derive an extinction map of the region, finding a median value of Ay = 4.33'" in 
the center of the association, decreasing toward the north-west. In the color-magnitude diagrams, the shape of 
the distribution of main sequence stars is compatible with the presence of an obscuring cloud in the foreground 
at ~ 850 + 25 pc from the Sun. 

Subject headings: stars: formation, stars: pre-main sequence, stars: color-magnitude diagrams (HR diagram), 
catalogs 



1. INTRODUCTION 

In recent years our knowledge of the star formation pro- 
cess and the early phase of stellar evolution has increased re- 
markably. This progress has been driven by a large number 
of deep spectroscopic and photometric observations at high 
spatial resolution of star forming regions in our Galaxy and 
in nearby galaxies. High performance telescopes, such as the 
NASA Hubble and Spitzer Space Telescopes, and more re- 
cently the ESA Herschel Space Observatory, allow us to study 
the formation of low- and high-mass stars, the evolution of cir- 
cumstellar disks, and the collapse of protostellar cores in un- 
precedented detail. X-ray observations with space telescopes 
such as Chandra and XMM-Newton have also been crucial, 
unveiling the populations of star forming regions down to sub- 
solar masses, and probing the high energy processes in young 
stars. 

Despite this progress, few single star forming regions in our 
Galaxy have provided the opportunity to study star and planet 
formation over a large range of stellar masses and in presence 
of very massive stars, which affect the star formation pro- 
cess and the evolution of their parental cloud mostly thank to 
their intense ionizing flux. The study of such massive young 
clusters, in fact, is hampered by their large distance from the 
Sun, with only a few moderately massive clusters closer than 
1.5 kpc. 

One exception that has only fairly recently been recog- 
nized is the Cygnus OB2 association ( Reddish et al.. .1967) . 
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in the Cygnus-X giant molecular cloud. It is the richest and 
most massive OB2 association within 2 kpc of the Sun. It 
harbors a large number of massive stars and an extensive 
young population of low-mass pre-main sequence stars, likely 
formed in different episodes. The first spectroscopic stu dy de- 
voted to the massive members of this region ( Reddis h et al.L 
119671) identified about 300 OB members, and in subsequent 
studies this estim ate has increased. In a 2MASS study, 
iKnodlsede^ ( l200tf ) used statistical arguments to suggest the 
presence of 2600 OB stars and 120 O stars and to estimate a 
total mass of (4 - 10) x IO'^Mq. This mass led the author 
to the conclusion that Cyg OB2 is a young globular clus- 
ter in the Milky Way. A near-infrared spectroscopic survey 
by IComeron et al.l (l2002h found slight ly fewer O stars but 
confirmed this co nclusion. While Ha nsonI (120031 see also 
iDrew et al] I2008 f) finds a somewhat lower total mass. The 
enormous scale of Cyg OB2 is exemplified by it con tainin g 



two of the few 03 stars known in ou r Galaxy (jWaibornl|1973h , 
together with B supergia nt stars (.Massev & Thompsonlll99ll : 
iNegueruela et al.L 120081). Other massi ve objects, such as 



Wolf-Rayet stars ('Niemela et al.', '1998") and the Luminous 
Blue Vai-iable G79.29-h0.46 (Higgs et al., 1994) are found in 
the Cygnus-X complex in the proximity of the association. 

The first studies of the intermediate- and low-mass popu- 
lation of the associa tion counte d several thousands of can - 
didate members tReddish et al.L 119671: iKnodlsedeii 120001) . 
but most recent and reliable estimates based on Chandra 
X-ray observations of the central cluster indicate a pop- 
ulation of 1000-1500 low -mass m embers down to subso- 
lar masses ( Albacete Colo mbo et al .. 2007; Wright & Drakel, 
i2009.). among which th ere are several strong Ha emitting ob- 
jects dVink etaUlIOOSi) . 

The first e stimate of the distance to Cyg OB2 was 2.1 kpc 
(iReddisheta l.. 1967). Massey & Thompson (1991) found 
1.8 kpc based on a combined photometric and spectroscopic 
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study, while the most recent estimate based on spectroscopy 
of the most massive members of Cyg OB2 is 1 .45 kpc 
( lHansonll2003h . iRygl et alj (1201 Ih reported a distance of the 
Cygnus-X complex of 1400 pc, derived from the trigonomet- 
ric parallaxes of five star forming region in the complex. We 
adopt a distance of 1450 pc in this work. Cyg OB2 then lies 
behind the nebulosity associated with the Great Cygnus Rift 
and, despite its relative proximity to the Sun, the association 
is significantly absorbed. 

A picture has emerged in the last few years of different 
episodes of star formation in and around the main center of 
Cyg OB2. iHansonl (flOOSh noted the presence of both mas- 
sive stars younger than ~ 2 Myr, and of B stars that appeared 
older for the adopted association distance, or would need to 
be in the foreground. iDrew et aTl 12008) identified a popula- 
tion of A-type stars south o f the m ain association with an age 
of 5 - 7 Myr. Since then, Wright et al. ( 2010) have shown 
that the NIR color-magnitude diagram of the region indicates 
both a ~ 3 Myr and a ~ 5 Myr ol d population of stars down to 
stellar masses of ~ 1 Mq, while IWright et al] (120 12h discov- 
ered a population of candidate proplyds^ southward Cyg OB2 
suggesting that new episodes of star formation occurred. 

In this paper we present the deepest and highest spatial res- 
olution optical observations of Cyg 0B2 to date. In sections 
|2] and [3] we describe the instrument setup, the observations, 
the data reduction, and the photometric calibration in the AB 
system; in section |4] we describe the final catalog, consist- 
ing of 64157 optical sources, among which 38300 have good 
photometry in all three bands. In section|5]we cross-correlate 
the OSIRIS catalog with the IPHAS data and with the newly 
released SDSS DR8 data of this region, and with the exist- 
ing Chandra/ AClS-l catalog of Cyg OB2. Finally, in sections 
|6] and [T] we present the color-color and color-magnitude di- 
agrams and we use them to study the extinction and age of 
association members and the distance of the foreground neb- 
ulosity. In a forthcoming paper, the optical catalog described 
in this paper will be combined with new deep X-ray obser- 
vations of Cyg OB2 (the 1 .08 Msec Chandra Cygnus OB2 
Legacy Survey, RI. J. Drake^), to obtain the deepest list of 
stars associated with this unique star forming region. 

2. INSTRUMENT SETUP AND OBSERVATIONS 

The observations were performed in the r', /', and z' filters 
with the Optical System for Imaging and low Resolution Inte- 
grated Spectroscopy (OSIRIS), mounted on the 10.4 m Gran 
Telescopio CANARIAS (GTC) of the Sp anish Obseryatori o 
del Roque de los Muchachos in La Palma (Cepa e t a^l2000^ . 
The OSIRIS detector comprises two Marconi CCDs, each 
with 2048 X 4096 pixels, separated along the long side by 
a 72 pixel gap. The total Field of View (FoV) is 7.8' x 8.5' 
(7.8' X 7.8' unvignetted), corresponding to a spatial scale of 
0.127" / pix. Such a small spatial scale is necessary to resolve 
very crowded regions such as Cygnus 0B2. 

OSIRIS provides different binning and frame readout 
speeds, each corresponding to different values of readout 
noise (RON) and detector gain^. Our observations have been 
performed with the 1x1 binning mode and with different 
readout speeds: lOOkHz/HN; 500kHz/HN and 500kHz/LN''. 

* Stars withjMxrtSpl^netaij disks surrounded by a photoevaporating enve- 
lope of gas (O'dell &l¥en, 1991) 

^ http://www.cygob2.org/ 

* on-line manual at http://www.gtc.iac.es/en/pages/instrumentation/osiris.pl 
' HN: High Noise, LN: Low Noise 



Table [T] summarizes the seeing conditions, configurations and 
observed Cyg OB2 and standard fields for each night of ob- 
servation. The seeing conditions were very good for the first, 
fourth and fifth nights; acceptable for the other nights. 

Fig. [Ushows a 1 .47° X 1 .47° IPHAS Ha image of Cyg OB2. 
The OSIRIS observations are arranged in a 5 x 5 mosaic in 
order to cover the 41' x 41' central part of the area imaged 
in X-rays for the Chandra Cygnus OB2 Legacy Survey (large 
box). Each OSIRIS field is labeled with a number ranging 
from 1 to 27. In the rest of this paper, we will refer to the 
fields using these numbers. Fields 16 and 17 are not shown in 
Fig.[T]since they were observed in very bad seeing conditions 
on 8/19/2009 and they have not been analyzed. Each field was 
observed three times per filter, for a total of 225 images, with 
small offsets in RA to cover both the gap between the CCDs 
and vignetted area. In this way, we obtained full coverage 
of the field in the RA direction, with some gap in the DEC 
direction due to CCD vignetting. The total exposure time in 
each field and in each filter was 80 sec, which allowed us to 
observe stars down to r' ~ 25'". The small inclined boxes 
in Fig. [H show the area previously observed wi th Chandra 
(iButt et al.Ll2006tlAlbacete Colombo etalll2007h . 

3. DATA REDUCTION AND PHOTOMETRY 

Overscan, bias, and flat field corrections were made using 
the IRAF CCDPROC task of the MSCRED package. No 
correction for dark current was been made since it does not 
significantly affect the OSIRIS detectors. 

Source detection and Point Spread Function (PSF) photom- 
etry were performed using the DAOPHOT/ALLSTAR pack- 
ages (lStetsodlT987llT99l . For source detection we used a 
threshold ranging from 3.5 to 4.5 times the standard deviation 
of the average sky emission level (ctsj.,,), and the stellar pro- 
files were fitted with a Gaussian function with a Full Width at 
Half Maximum (FWHM) ranging from 0.9" to 1.5", depend- 
ing on the seeing conditions. The PSF photometry has been 
performed fitting the observed stellar profiles with different 
PSF models, described by either a Gaussian, Moffat, Lorentz, 
or Penny functions. We also accounted for a possible linear 
variation of the PSF model with the plate coordinates. 

In order to combine the different observations of each 
star, we calculated the plate transformations between the 
different images of the same field, by using DAOMAS- 
TER/DAOMATCH procedures dStetsonl 119931) . adopting a 
linear translation of the coordinate systems. Self-consistent 
sets of PSF stellar magnitudes and positions between all the 
images of the same field have been obtained using ALL- 
FRAME ([stetson, 1994) and the plate transformations calcu- 
lated wifli DAOMASTER/DAOMATCH. We calculated the 
aperture correction to the photometry using the DAOGROW 
routme dStetsonl [1990), which uses the growth - curves 
method, and a detailed model for the stellar profile. 

3.1. Photometry and calibration 

Photometric calibration in the SDSS riz photometric sys- 
tem was carried out using the CCDSTD/CCDAVE/TRIAL 
procedures ( Stetson, 2005). These three routines calculate the 
coefficients of the chosen photometric solution, then apply the 
transformation first to the standard stars together with a se- 
lected sample of target stars, and then to the whole catalog. 
We calibrated the catalog in the SDSS riz photometric system 
since the calibrated magnitudes of the observed standard stars 
were taken from the eighth data release of the Sloan Digital 
Sky Survey (SDSS DR8. .Aihara et al.„.2011,) . 
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Table [T] summarizes the observations of the standard fields. 
Unfortunately, SDSS DR8 data are not available for the field 
GD71, meaning that the standard stars necessary for the cali- 
bration of the images taken on August 18th were not available. 

We calibrated the observations of Cyg OB2 taken on the 
other nights (all the fields with the exception of 6 and 7) using 
the SDSS photometry of the standard stars in the observed 
standard fields. We used a standard photometric solution: 

Ok^ Mk+ZPk+AlkXCT +A2kXQ (1) 

Here, k is one of the photometric bands (r, /, or z). Ok is 
the instrumental magnitude in the k band, Mk is the calibrated 
magnitude, ZP^ is the photometric zero-point, CT is a first or- 
der color term, and Q is the airmass. Since the standard stars 
have not been observed over a wide range of airmass, we set 
A2 equal to the extinction in r' measured during the observing 
nights by the Carlsberg Meridian Telescope'*', converted in /' 
and z' using the extinction curve derived by King 1985". We 
calculated A 1 from the stars observed in both CCDs, and then 
two diff'erent zero-points from the CCDs separately. Table |2] 
shows how the coefficients used for the photometric calibra- 
tion vary in the 5 observing nights. Residual zero-point and 
color corrections have been applied for e ach field after match- 
ing the GTC and SDSS DR8 (Sect. l5H l. 

Fields 6 and 7 have been calibrated by taking advantage of 
the fact that field 8 has been observed both on August 18th 
and 20th. We simply selected a sample of good stars compat- 
ible with the MS falling in this field, and then we used their 
calibrated photometry to propagate the photometric solution 
in fields 6 and 7. 

The final photometric solutions reproduce properly the cal- 
ibrated magnitudes of standard stars from the instrumental 
magnitudes, for the whole range of stellar magnitudes, col- 
ors, positions and airmass. In all cases, the distributions of 
the residuals (the difference between the calibrated and stan- 
dard magnitudes for the standard stars) have medians close to 
zero and RMS deviations of about, or smaller than, 0.035. 

The astrometric solution was found using the IRAF tools 
CCXYMATCH, SETWCS SiudSKYPIX, adopting the 2 Mi- 
cron All Sky Survey (2MASS) Point Source Catalog (PSC, 
ICutri et al.. 2003) as astrometric reference and with a tanx 
projection of the plate coordinates onto the celestial system. 

4. THE FINAL CATALOG 

Table|3]shows part of the final OSIRIS catalog of Cyg OB2. 
For each star it provides an identification number (column 1), 
the celestial coordinates (columns 2 and 3), magnitudes and 
errors in r, /, and z (columns 4 to 9), and the and SHARP 
parameters. The former is the chi-square derived in the PSF 
fitting phase, the latter is a para meter describing the shape of 
the stellar profile (IStetsonl|l987 ). These parameters lose their 
full statistical meaning when the PSF models from different 
images of the same star are combined in the ALLFRAME 
phase, but they are still useful to select stars with good and 
bad photometry. 

The merged catalog consists of 64157 sources. As ex- 
plained in Table [3] stars with and SHARP parameters equal 
to have been detected only in one band (z' usually). The 
magnitudes of these stars are not calibrated properly (they do 
not have an available color). Since they might be useful for 

http://www.ast.cam.ac.uk/ioa/research/cmt/camc_extinction.html 
' ' http://www.ing.iac.es/Astronomy/observing/manuaIs/ps/tech_notes/tn03 1 .] 



future multi-wavelength studies involving the OSIRIS catalog 
they are retained in order to preserve the locations of known 
faint optical sources. These stars can be easily selected by 
requiring - (there are 50 of these null - stars). 

Stars with bad PSF fitting or with non-circular profiles can 
be easily selected by requiring > 5 or \sharp\ > 8. A total 
of 1815 stars (2.8% of the whole catalog) satisfy these crite- 
ria. Most of them are sources close to the CCD edges and 
corners, diffraction spikes created by bright stars, or saturated 
sources with irregular profiles. The total number of sources 
with both good shape and but with errors in colors larger 
than 0.15'" is 23992, or 37.4% of the whole catalog. Finally, 
the sources with good photometry in all three bands number 
38300, or 59.7% of the catalog. We will call these sources 
good - photometry sources hereafter. The saturation limits 
have been estimated, by inspection of the stellar profiles, to 
be r' = 16"', ;' = 15.5"', and z' = 15'". Fig.|2]shows the spa- 
tial distribution of the good - photometry stars of the OSIRIS 
catalog. The gaps in declination are due to the fact that stars 
falling in vignetted areas do not have the properties required 
for the good- photometry credential. The two Chandra ACIS- 
I pointings of Cyg OB2 presented bv lAlbacete Colombo et alj 
(2007) and Butt et al. (2006) are also shown. In Fig. |2] there 
is an evident increase in the density of the good - photometry 
stars from the south to the north, indicating a variation of the 
visual extinction in the south-north direction. In fact, the com- 
parison of Fig.|2]with infrared images of Cyg OB2 (not shown 
here) reveals that the area to the north-west, which is that 
with the highest visible stellar density, corresponds to a cavity 
in the molecular cloud where a large number of background 
sources can be observed (Red dish et aUlT967h . 

Fig. [3] shows the photometric errors vs. the magnitude in 
the three bands. In each panel, the horizontal cr - 0.1"' line 
shows the depth of the photometry of the good - photometry 
stars, reaching r' ~ 25"', /' ~ 24'", and z' ~ 22.5"'. 

5. CATALOG CROSS -MATCHING 

5.1. Cross-correlation with the SDSS DR8 catalog 

The data from the Eighth SDSS Data Release cover the en- 
tire OSIRIS field, with the exception of the central area corre- 
sponding to the location of the most massive stars. The SDSS 
catalog is therefore not useful to study the central cluster, and 
it is also 3 magnitudes shallower than the OSIRIS catalog, but 
it does provides useful information such as the magnitudes 
in u and g bands important for studying stellar photospheric 
properties, and important means for verifying the photometric 
calibration. 

We matched the OSIRIS and SDSS DR8 data in two steps. 
The first step was aimed to define the correct matching radius, 
by experimenting with increasing radii from 0.1" to 2" and 
comparing the differential distributions of observed, real, and 
spurious matches at increasing matching radii. The number of 
spurio us coincidences has been calculated as in lDamiani et alJ 
(120061) . taking into account that the two catalogs are corre- 
lated; the total of real matches is given by the difference be- 
tween the observed and spurious coincidences. We adopted 
a final matching radius of 0.3", which is the highest value 
at which the differential distribution of spurious matches is 
negligible (few percent) compared with that of real matches. 
At this radius, we found 10790 matches between the OSIRIS 
and IPHAS catalogs, with 54 expected spurious matches and 
8 multiple matches where a source of one catalog is matched 
with more than one of the other. After correcting for residual 
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zero-points and color dependence in each OSIRIS field, we 
achieved negligible median offsets between the OSIRIS and 
SDSS magnitudes smaller than 0.003"'. 

5.2. Cross-correlation with the IPHAS catalog 

Another catalog in three optical bands (r', /', and Ha) of 
the Cyg OB2 region can be extracted from the observations 
taken with the Wield Field Camera (WFC) on the 2.5m Isaac 
Newton Telescope (INT) f or the INT Photometric Ha Survey 
(IPHAS. iDrew et al.ll200l . A total of 11995 IPHAS sources 
fall in the OSIRIS field. Our new OSIRIS data are about 5 
magnitudes deeper in r and / than the IPHAS catalog. How- 
ever, the latter is 4 magnitudes brighter than our saturation 
limit in the r band. A combination of both catalogs allows 
us the study of stars over a wide range of magnitude, from 
r - 12'" to r = 25"'. One complication is that, as described 
in Sect. 13.11 our OSIRIS catalog has been calibrated in the 
SDSS photometric system (AB), while the IPHAS catalog is 
calibrated in the Vega system. Suitable color transformations 
between the two catalogs have to be found in order to convert 
the IPHAS colors into the AB system. We adopted an empir- 
ical approach whereby the color transformation was obtained 
from sources common to both catalogs. 

To match the IPHAS data with the SDSS-hOSIRIS cata- 
log, we adopted the same procedure described in the previ- 
ous section, using a matching radius equal to 0.5", and find- 
ing 111 99 matches between the OSIRIS and IPHAS catalogs 
(17 expected spurious matches), divided in 6104 sources with 
OSIRIS, IPHAS and SDSS counterparts and 5095 with only 
OSIRIS and IPHAS counterparts. 

A set of photometric transformations between the SDSS 
and WFC photometric system s has been already found by 
iGonzalez-Solares et al.1 (1201 ll) : 

r'wFC ^rsDss -0.144 + 0.006 x(g- r)soss (2) 
i'wFC - isDss - 0.41 1 - 0.073 x (r - i)sDss 0) 

with a weak color term in r band given by the similarity 
of the two used r filters. Since the SDSSh-IPHAS data of the 
Cygnus field c over a larger range of colors ( g - r < 4"') than 
those used bv iGonzalez-Solares et al.l (1201 11) (g - r < 1.8™), 
we reviewed these transformations, selecting stars with errors 
in the relevant colors smaller than 0.15'" and in the magnitude 
offsets smaller than 0.03'". We also removed in four iterations 
the sources with an offset more deviant than lia-offset from 
the value expected from the linear interpolation of the data. 
The result ing transformations are slightly different than those 
found bv lGonzalez-Solares et al.1 (1201 Ih : 

r'wFC = rsDss - 0.105 - 0.034 x{g- r)soss (4) 
i'wFC - isDss - 0.351 - 0.071 x (r - i)sDss (5) 

mostly for the presence of a significant color term in the 
r band, which is very likely induced by the small difference 
in extinction at th e central wavelen gths of the two r filters 
(see, for instance, lO'DonnelH (119941) ). whose effect is domi- 
nant only for very red objects. 

To convert the WFC magnitudes into the SDSS system we 
used the r - i color in both magnitudes, expecting to find 
a significant color dependence like those in Eq. U and |5] 
Fig. |4] shows the magnitude offsets of rgsmis - i'wFC 
iosiRis ~i'wFc ^ function of the OSIRIS r-i color. The stars 



used in this diagram are the good - photometry sources cor- 
responding to single matches between the two catalogs, with 
r and r' magnitude in the 14'" - 18'" range, and with the error 
of the difference between the OSIRIS and IPHAS magnitudes 
smaller than 0.03"' (which is the median difference error for 
all the matched sources). We also rejected in four iterations 
the stars whose offsets are more deviant than 'ia-^nffgi-ence from 
the value expected from the linear fit to the data (the fit was 
performed for each iteration). The coefficients obtained from 
the linear interpolation in the last iteration are shown in the 
top of each panel in Fig.|4] The transformations between the 
OSIRIS photometric system (the AB system) and the IPHAS 
one (the Vega system) that we obtained are then: 

rosmis - r'wFc = 0.072 + 0.085 x (r - i)osiRis (6) 
iosiRis - i'wFC = 0.331 + 0.072 x (r - Oos/r/s (7) 

where rgsmis and iosms are the OSIRIS and r^^^ and 
i\^P(^ the IPHAS magnitudes. Eq.|7]gives the following color 
transformation: 

('■' - HwFc = 0.260 + 0.987 x (/ - /')os/r« (8) 

To verify these transformations, we used them to con- 
vert the expected r' - i' color of a G2 V star. We adopted 
as the color value in the Vega system dDrew et aTl 12005) . 
(r' - i')vp„n = 0.384"', and in the AB system the value found 
bv iFukugita etaH (l20Toh . (r'-i%B = 0.127"', equal to 

~ Oab - 0.126™ by using the conversion suggested in the 
SDSS DR8 website' . Using Eq. [8] to convert the color of 
a G2V star from the AB system to the Vega system, we find 
0.384"', reproducing properly the Vega r - i color of a solar 
type star. 

5.3. Cross-correlation with the Chandra catalog 

Thanks to the high level of X-ray emissio n characteriz- 
ing young PMS stars (e.g. iMontmerlH 11996b X-ray selec- 
tion provides a powerful means to isolate the cluster pop- 
ulation from the myriad foreground and background stars. 
The combination of our optical catalog with X-ray data al- 
lows us to identify and study the cluster loci in the optical 
diagrams. Two existing Chandra ACIS-I observations par- 
tially cover the field observed with OSIRIS. The first one (P.I.: 
Flaccomio) centered on Cyg OB2 (ff = 20 : 33 : 11.0 and 
6 - +41 : 15 : 10.00) with a total exposure tinieof 99.7 A;i 
was first analyzed by Albacete Colombo et all (l2007h . The 
second one (PL: Butt) covers a field to the north-west, cen- 
tered on Of = 20 : 32 : 07.0 md 6 ^ +41 : 30 : 30.00, with 
an exposur e time of 49.4fc^ . Both observations have been an- 
alyzed by Wright & Drakd (120091) . who derived a catalog of 
1696 X-ray sources complete down to 1 Mq. The FoVs of the 
two Chandra observations are shown in Fig. [T] 

We first performed a preliminary match between the optical 
and X-ray catalogs by using a 2" matching radius, in order 
to determine and correct systematic position offsets. Since 
the PSF in the ACIS-I detector (and the resulting instrument 
sensitivity) degrades from the center of the field outward, for 
the final match we used a variab le ma tching radius. Moreover, 
the procedure adopted in Sect. l5.2l to derive the distribution 
of spurious matches at different matching radii requires the 
matched catalogs to be uniformly distributed in the common 

http://www.sdss.Org/dr7/algorithms/jeg_photometric_eq_drl.html#usno2SDSS 
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area. This is a good approximation for the optical catalogs, 
but, owing to the variable sensitivity of the Chandra FoV, the 
spatial distribution of Chandra sources can not be considered 
uniform. 

To account for the varying PSF, each X-ray source has 
been matched with close optical sources by using a match- 
ing radius proportional to its positional uncertainty: r,„a,ch - 
A X o-pos- The proportional constant A is common for every 
X-ray source, and it has been found by comparing the distri- 
butions of spurious and real matches obtained with 20 values 
of A ranging from 0.1 to 5. For each value of A, we per- 
formed 4 matches, applying a rigid translation to the X-ray 
catalog in order to derive the number of expected spurious 
matches Nspmious {A). The total number of matches Ntotai {A) 
has been derived by matching the X-ray catalog without trans- 
lations, and the real number of matches Nreai {A) - N,oiai {A) - 
N spurious {A). Comparing again the distributions of real and 
spurious matches, we set the value of A equal to 1.5, this be- 
ing the largest value at which the differential distributions of 
spurious coincidences is negligible with respect to that of real 
matches. With this definition of the individual matching ra- 
dius, we obtained a total of 1407 matches, with 22 expected 
spurious matches. This result is a great improvement over the 
catalog published by Wright & Drakd (2009), who matched 
the X-ray and existing IPHAS catalog of this region finding 
750 coincident sources. 

6. VISUAL EXTINCTION IN CYG OB2 
6.1. Extinction in the direction of the association 

Despite the proximity of Cyg OB2 compared with other 
massive clusters, it is affected by a large visual extinction due 
both to the dust associated with the Great Cygnus Rift in the 
foreground and its parental cloud. The most recent estimates 
of the visual extinction in the direction of Cyg OB2 are those 
of lSale et al.l (120091) . who studied how A v varies with distance 
in the direction of the association and found an increase of Ay 
from 2™ to 5"* in the region between 1 kpc and 2 kpc away; 
iDrew et aTl (120081) . who found that the cluster locus in the 
r' - Ha vs. r' - i' exhibits an extinction between Ay = 2.5™ 
and Ay - 7™, with a larger extinction in the area surrounding 
the association; and Wright et al. (2010) . who found a median 
Ay = 7.5™ in the central cluster and Ay - 5.5"' in the north- 
west field observed by Chandra. 

Since colors are not affected by distance, the color-color r-i 
vs. / - z diagram is suitable for estimating the visual extinc- 
tion affectin g Cyg OB2 members, by fitting the isochrones 
from Siess et al I (l2000t) to the cluster locus. However, the 
photometric data of the isochrones are in the Johnson-Cousin 
UBVRI photometric system, so they must be converted us- 
ing proper transformations from the U BVRI to the ugriz sys- 
tem. Several transformation have been published. In order to 
choose the most reliable for our case , we compared the r - i 
and / - z colors of the ZAMS from ISiess et all (120001) con- 
verted using different transf ormations w i th thos e of the main 
sequence stars observed by ICovev et"aLl (l2007h . as shown in 
Fig. |5] The transformation that rep roduces better the col - 
ors of the Covey MS stars is that from Fukugita et all (Il996h . 
This transformation converts the f/Z?y/?/ photometry into the 
u'g'r'i'z' system, so it is necessary to apply also the trans- 
formation between t he u'g'r'i'z' and th e ugriz systems. The 
transformation from lSmith et al.l ( 120021) fails to re produce the 
colors of late stars. The transformations from iJordi et al] 
(12006.) and .Rodgers et al.. (2006.) transform the Siess ZAMS 



into a straight line. The former is a better appr oximation of 
the irr egular shape of the Main Sequence from ICovev et"an 
( 120071) . but the difference between the two transformat ions is 
small. For these reasons, we mainly adopt the Fukugit a et al.l 
('1996) transformations. 

Fig. |7] shows the r - i vs. / - z diagram of all the 
good - photometry stars in the OSIRIS FoV. Since almost all 
the optical sources with X-ray counterparts are young clus- 
ter members, they define a well determined locus in the di- 
agram. In Fig. Ill we also show the colors predicted by the 
3.5 Myr isochrone (ISiess e t al., 2000), which isjhejnean age 
of the stars in the center of Cyg 0B2 ( Wri ght et al.l 120101) . 
transformed into the ugriz system and redde ned by using the 
extinction law derived by O'Donnelll ( 119941) . From the com- 
parison of the observed colors with this isochrone we find 
the extinction of the cluster members to range mainly from 
Ay = 2™ to Ay = 6™. A few X-ray sources show higher ex- 
tinctions, being largely embedded cluster members or back- 
ground sources. A large group of sources lies to the left of 
the 3.5 Myr isochrone with Ay = 2"', and is separated from 
the cluster locus by a small gap with a smaller density of stars 
compared with the rest of the diagram, being mostly in the 
foreground. 

In principle, the cloud responsible for the extinction affect- 
ing Cyg OB2 can be in the foreground or associated with the 
cluster However, there are several indicatio ns that Cyg OB2 
is not yet embedded in the parental cloud ( Schneider et al.L 
l2006h and it lies behind the Cygnus Rift, which is responsible 
for a steep rise of the visual extinction along this direction. 
Since the optical colors at short wavelength are largely af- 
fected by extinction, it is possible to estimate the distance of 
the nebulosity which is mainly responsible for the rise of vi- 
sual extinction from the color-magnitude diagram g vs. g - i. 
The left panel of Fig. |6] shows 100 ZAMS from isiess et alJ 
(2000) in this color-magnitude diagram, each drawn adopting 
a visual extinction Ay = 0.5 x distikpc), typical of lines of 
sight crossing an uniform and loose distribution of interstellar 
medium, and distance ranging from 500 pc to 5500 pc. As ex- 
pected, increasing distance and extinction the ZAMS populate 
the lower part of the diagram. The right panel of Fig. |6]shows 
the observed gMS. g- i diagram of the OSIRIS-SDSS sources 
with good photometry (small dots), with those detected in X- 
rays marked by l arge d ots. The dashed lines are isochrones 
from ISiess et al.l (12000 1) drawn with the distance adopted for 
Cyg OB2 (145 pc) and the average extinction of Ay = 4.2"' 
(see Sect. 16.21 ). Comparing directly the two panels, it is evi- 
dent that the observed diagram is less populated by faint blue 
sources than what expected from the left panel. This is due to 
the presence of the nebulosity along the line of sight, which 
causes a steep rise of the reddening of the stars behind it. The 
distance of this cloud can be estimated, then, by finding what 
are the values of Ay and distance at which the ZAMS fits the 
blue end of the observed distribution of stars in the g vs. ^ - / 
diagram. The value of distance at which we obtained the best 
fit is 850 ± 25 pc, with Ay - 1™. This value is on ly slightly 
smaller than the distance found bv^ Sale et al.l (12009 ). who es- 
timate that the extinction along this line of sight starts to in- 
crease from Ay = 2™ (~ 1 kpc) to Ay = 5™ (~ 2kpc). Our 
estimate also confirms that the nebulosity responsible for the 
high visual extinction is not associated with Cyg OB2, but is 
in the foreground. 

6.2. Individual extinction of candidate cluster members 



6 



Taking advantage of the fact that the 3.5 Myr isochrone in 
the r - i vs. / - z diagram can be transformed int o a straight 
line with the transformations of iJordi et alj (120061) with good 
approximation, it is possible to calculate the individual ex- 
tinction of the candidate cluster members from the the dis- 
placement of the X-ray sources from the 3.5 Myr isochrone 
drawn with Ay - 0"' along the reddening vector. Indeed, we 
verified that contributions to these displacements from the ex- 
pected range of ages of Cyg OB2 cluster members are negligi- 
ble. In this calculation, we excluded the few sources more red 
than the colors of a 3.5 Myr old B2.5 star at different extinc- 
tions. The variation in apparent extinction could be caused by 
a patchy distribution of gas and dust along the line of sight, 
or by intrinsic circumstellar gas such as a protoplanetary disk. 
The presence of a circumstellar disk is only expected to sig- 
nificantly affect the optical colors of stars when intense accre- 
tion is ongoing, or for a disk inclination larger than 80° with 
respect to the line of sight (Guarcello et al., 2010). For most 
sources then it is likely that the variations in extinction corre- 
spond to an inhomogeneous distribution of gas and dust in the 
cluster environment. 

Fig.[8]shows the distributions of individual extinctions cal- 
culated for the optical sources with X-ray counterparts falling 
in both Chandra fields (upper panel) and in each field (lower 
panels). We have not considered sources with Ay < 2"\ since 
they are more likely foreground stars (a total of 56 objects). 
Sources in the southern Chandra field (centered on Cyg OB2) 
have a median extinction equal to 4.33"'. This m edian extinc- 
tion is smaller than the value found by Wright et al.l (1201 Oh 
with an isochrone-fitting technique, but it is withi n the range 
of ext inction between 1 kpc and 2kpc fou nd by ISale et afl 
(l2009h (2 < Ay < 5) and lDrew et al! (l2008h (2.5 < Ay < 7). 
In the histogram (middle panel) the large differential redden- 
ing affecting the central cluster is also evident, since the dis- 
tribution does not show an evident peak and the bins from 
Ay = 3.0"' to Ay - 5.0"' are almost equally populated. The 
extinction distribution in the north-west field is quite different, 
with a median value equal to 3.21"", and a well defined peak 
with no long tail toward larger extinction. This indicates that 
the density of the obscuring material changes significantly in 
the direction of the two fields. 

To better visualize how the extinctions vary across the field, 
the left panel of Fig. |9] shows a gray-scale map of the visual 
extinction, where the Ay value characterizing each bin is the 
median value of the extinction of the sources falling in the bin. 
The median number of sources in each bin is 1 1, with the less 
populated bins lying along the border of each ACIS field. The 
locations of the known O stars are also indicated. In the north- 
west field the extinction distribution is almost uniform, while 
in the southern field the extinction increases in the center-west 
direction. 

The right panel of Fig. |9] shows the density map of the can- 
didate cluster members (optical-i-X-ray sources with Ay > 
2'"). In the southern field the density of cluster members, 
both low-mass and O stars, clearly peaks in the center of the 
field, corresponding to one of the most extinguished regions. 
This could hint that, despite the large content of massive stars, 
Cyg 0B2 is still embedded in the parental cloud. However, 
the CO map shown in Schneider et al. (2006) suggests that the 
cluster has created a cavity in the parental cloud, as is com- 
monly observed in massive clusters in our Galaxy with an age 
of a few Myrs. On the other hand, in this CO map small re- 
gions of higher density lie close to the positions of the two 
groups of massive stars and could be responsible for the peak 



of extinction evident in Fig.|9] The selection of the members 
with disks, with the identification of those members still em- 
bedded in dense regions of the parental cloud, will shed some 
light on the possibility that embedded sources and dense intr- 
acluster medium is still present in the center of Cyg OB2. 

7. COLOR-MAGNITUDE DIAGRAMS OF CYG OB2 
MEMBERS 

Fig. [To] shows the three color-magnitude diagrams of the 
good - photometry sources in the whole OSIRIS field. Stars 
with X-ray counterparts are highlighted. In all the diagrams, 
the X-ray sources define a well delimited cluster locus. By 
drawing the isochrones younger than 10 Myr with an extinc- 
tion of Ay = 4.2"', the cluster loci in all the color-magnitude 
diagrams mostly lie between the 5 Myr and 1 Myr isochrones. 
The ZAMS, drawn adopting the distance of the foreground 
cloud of 850 pc and an extinction of Ay = 1"', is well fit- 
ted to the blue end of the locus of main sequence foreground 
stars, indicating that these parameters are appropriate. Most 
of the X-ray sources bluer than the 10 Myr isochrone are those 
with extinction smaller than 2"' that are very likely in the fore- 
ground. 

The r vs. r - i diagram in Fig. [TO] is completed using 
the IPHAS data, converted into the SDSS photome tric sys- 
tem by using the transformation we found in Sect. 15.21 all 
the three diagrams are also completed with the SDSS data. 
Note that, thanks to the depth of the OSIRIS observations, we 
now have optical photometry for a significant number of faint 
low-mass candidate cluster members: 484 X-ray sources with 
r > 21.5"' out of a total number of 1250 sources with both 
X-ray and optical counterparts with good photometry. This 
will be very important for future studies based on the selec- 
tion of low-mass cluster members and the effects of massive 
stars on their early evolution. The masses listed on the right 
of the diagram, corresponding to the 3.5 Myr isochrone (rep- 
resenta tive of the median age of the central cluster as noted 
earlier; Wright et al.L 12010 ) demonstrate that the grasp of the 
OSIRIS photometry is deep enough to observe stars with sub- 
solar mass at the distance of Cyg OB2. 

Among the IPHAS sources with X-ray counterparts, those 
with r - i ~ 0.5"' and 12"' < r < 16"' are compatible with B 
stars with only moderate extinction. Most of them are in the 
north-west field, where the extinction is lower in general, as 
might be expected. A group of X-ray sources withr-/ ~ 0.1"' 
and 13.5"' < r < 12"' are compatible with being foreground 
objects. This is also consistent with their positions in the r - i 
vs. i - z and r' - Ha vs. r' - i' diagrams and with their low 
neutral hydrogen column densities, Nh, deduced from their 
X-ray spectra that are compatible with extinctions Ay < 0.2"' 
d Wright et al.L 12010). Surprisingly, their spatial distribution 
is not sparse as expected for uncorrected foreground objects, 
but instead they all lie in the north-west field. 

8. CONCLUSIONS 

We have analyzed new GTC/OSIRIS optical observations 
in r'i'z' bands of a field of size 41' x 41' approximately cen- 
tered on Cyg OB2. The resulting catalog contains the photom- 
etry of 64157 optical sources, among which 38300 have good 
photometry in all three bands. The catalog reaches r' - 25"', 
5 magnitudes deeper than the existing IPHAS catalog and 3 
than the SDSS DR8 data. This Hmit corresponds to a 3.5 Myr 
star with M - 0.15 Mq at the adopted distance and extinction 
of Cyg OB2. The stellar density of stars with good photom- 
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etry varies significantly from north southward, with the max- 
imum density corresponding to a cavity in the cloud in the 
north-west part of the observed area. 

We cross-correlated the OSIRIS catalog with an existing 
Chandra ACIS-I catalog, of the central and north-west area 
of the field, in order to define the cluster locus in the dia- 
grams. This yielded 1407 optical sources with X-ray counter- 
parts, almost twice the number of X-ray sources with optical 
counterparts identified in previous works in these fields. Us- 
ing the stars with detections in both the OSIRIS and IPHAS 
catalogs, we derived a suitable color transformation between 
the IPHAS and OSIRIS photometric systems (Vega and AB, 
respectively). This transformation converts properly the col- 
ors of a G2V star from the AB to the Vega system. 

The r - i vs. i - z color-color diagram shows a clear cluster 
sequence which can be encompassed by a 3.5 Myr isochrone 
with extinction between 2"" and 6"". This extinction turns out 
to be the main range of extinction of cluster members. We 
also found that the mean extinction decreases from the central 
cluster northward: in the central field we found a median ex- 
tinction Ay =4.33™, with a large diff'erential reddening across 
the field; while in the north-west Chandra field we found a 
median extinction Ay = 3.21'", with a more peaked extinc- 
tion distribution. We derived an extinction map of the region 
that reveals evidence that the most obscured regions are also 
those with the largest density of low-mass and massive mem- 
ber stars. 



The color-color diagram also exhibits a gap between the 
foreground objects and the cluster locus that is induced by 
a steep increase of the visual extinction. This is due to the 
presence of a dense nebulosity along this line of sight. The fit 
of the zero age main sequence with the foreground population 
in the g vs. g - i diagram suggests a distance of 850 ± 25 pc 
for this cloud, which turns out to be in the foreground and not 
associated with Cyg OB2. 

This article is based on observations made with the Gran 
Telescopio CANARIAS (GTC), installed in the Spanish Ob- 
servatorio del Roque de los Muchachos of the Instituto de As- 
trofisica de CANARIAS, in the island of La Palma, and it 
also makes use of data obtained as part of the INT Photomet- 
ric Halpha Survey of the Northern Galactic Plane (IPHAS) 
carried out at the Isaac Newton Telescope (INT) (all IPHAS 
data are processed by the Cambridge Astronomical Survey 
Unit, at the Institute of Astronomy in Cambridge) and from 
archival Chandra/ACIS-I observations. MGG and NJW were 
supported by Chandra Grant GO0-11040X. JJD, VLK, and 
TA were supported by NASA contract NAS8-39073 to the 
Chandra X-ray Center (CXC) and thank the Director, Harvey 
Tananbaum, and the CXC science staff for advice and sup- 
port. DGA acknowledge support from the Spanish MICINN 
through grant AYA2008-02038. 
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TABLE 1 

Loo OF THE OBSERVAnONS 



Night 


Fields 


Standard Fields 


Configuration 


Seeing {r'J',z') (arcsec) 


8/1 1/2009 


11-15 


PG1545; SP0346 


200/LN 


0.90,0.94, 1.00 


8/18/2009 


6-8 


GD71 


500/HN 


1.29, 1.21, 1.28 


8/20/2009 


1-5 


PG0231; SA92 


500/LN 


1.29, 1.21, 1.28 


8/21/2009 


8-10 


PG2317 


500/HN 


1.18, 1.12, 1.15 


8/21/2009 


18-22 


PG2317 


500/LN 


1.08.0.99, 1.04 


8/22/2009 


23-27 


SA95 


50()/LN 


1.47, 1.43, 1.48 



TABLE 2 

Coefficients of the photometric caubration 



r 


i 


Z 


Left CCD zero-points 


3.58-4.00 


3.18-3.61 


2.76-3.22 


Right CCD zero-points 


3.55-4.02 


3.15-3.61 


2.80-3.29 


Color terms 


0.266-0.318 


0.065-0.087 


0.013-0.029 




J 



1- 


2 


1 

3\ 


4 


5\ 

- 1 


^6 

■ 


-A- 


• •■■tJ 


M— 

:s 


-It)' 



XL 



'i2 




13 



15 







-■20 


'2l 


22 


23 


'24' 


25 


26 


27- 



OSIRIS fields 



34:00.0 33:00.O 32:00.0 




31:00.0 20:30:00.{ 
. I . . . I , • 



Fio. 1. — IPHAS Ha image of a 1.47° x 1.47° area centered on Cyg OB2. The numbered black boxes are the fields observed with OSIRIS. The large inclined 
dashed box represents the area observed in the Chandra Cygnus 0B2 Legacy Survey. The small inclined boxes are earlier Chandra pointings in this region. 



TABLE 3 
Part of the electronic catalog. 



ID 


RA(J2000) 


DEC(J2000) 


r' 


(Tri 


/' 


cr,' 


z' 


o"z' 


>^ 


SHARP 


171 


308 6602478 


41.4834251 






24.193 




22.012 




0.00 


0.00 


172 


308 6459961 


41.4840431 


22.673 


0.012 


20.673 


0.005 


18.940 


0.006 


1.26 


0.50 


173 


308 6399231 


41 4851074 


24.659 


0.159 


22.822 


0.037 


21.072 


0.028 


1.19 


1.13 


174 


308 6605835 


41 4853439 


25.513 


0.327 


23.922 


0.103 


21.748 


0.126 


1.25 


4.56 


175 


308 6415100 


41 4851723 


22.106 


0.010 


20 688 


0.007 


19.425 


0.014 


1.16 


0.37 


176 


308.6547546 


41.4855270 






23.922 




22.405 




0.00 


0.00 


177 


308.6530151 


41.4855690 


21.426 


0.007 


19.808 


0.005 


18.920 


0.006 


1.20 


0.15 


178 


308.6368103 


41.4860077 


24.532 


0.122 


22.391 


0.021 


19.908 


0.009 


1.29 


1.32 


179 


308.6352539 


41.4861221 


24.496 


0.109 


22.598 


0.039 


21.137 


0.040 


1.21 


1.47 


180 


308 6576538 


41 4866714 


21.586 


0.007 


20.192 


0.005 


18.933 


0.011 


1.17 


0.10 


iol 


JVO.O'+l'+lOJ 






U.UjD 




n ni 5 




U.Ul J 


1 1Q 

1.17 


ft 1 9. 
U. 10 


182 


308.6518250 


41.4867249 


24.573 


0.087 


22.790 


0.057 


21.194 


0.036 


1.18 


0.46 


183 


308.6606445 


41.4868431 


25.838 


0.391 


23.777 


0.108 


21.720 


0.036 


1.21 


1.36 


184 


308.6435242 


41.4867668 


22.842 


0.017 


21.262 


0.009 


19.937 


0.014 


1.17 


0.42 


185 


308.6415710 


41.4869843 


24.671 


0.070 


22.630 


0.038 


21.058 


0.019 


1.15 


0.39 


186 


308.6609802 


41.4873161 


25.137 


0.289 


24.342 


0.213 


22.209 


0.168 


1.21 


1.32 


187 


308.6427612 


41.4873199 


23.695 


0.033 


21.560 


0.015 


19.805 


0.011 


1.19 


0.18 


188 


308.6546021 


41.4875298 


18.918 


0.005 


17.131 


0.001 


15.452 


0.001 


1.62 


-0.02 


189 


308.6382141 


41.4877434 






24.594 




22.874 




0.00 


0.00 


190 


308.6494446 


41.4878654 


24.681 


0.189 


24.275 


0.215 


21.945 


0.054 


1.18 


2.64 



If both;;^ and SHARP are equal to then the source was observed in only one optical band, 
and it was not possible to caUbrate its magnitude. 
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Fig. 3. — Photometric errors vs. magnitudes in r', i' , and z' for all the stars in the OSIRIS catalog. The horizontal line marks the cr = 0.1"' limit. 
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0,07a(± 0.002)+ 0.085(± 0,003)x(r'-i')„, 
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r'-i' (OSIRIS) 

Fig. 4. — Magnitude differences in rosiRis ~'''iphas ^"'^ 'OS iris ~''iphas ^^^sus the OSIRIS r-i color for the stars in common between the OSIRIS and IPHAS 
catalogs selected as described in the text. The continuous lines, whose coefficients are shown in the top of each plot, are obtained from a linear fit to the data. The 
dotted line in the top panel marks the median r offset. 




Fig. 5. — - r — i and i — z c olors of the Main Sequence stars with different spectral types observed ICovev et alj 420071) (dotted line), and those predicted by the 
ZAMS of lSiess et aP 120001 with different transformation from the UBVRI to the ugriz photometric system applied. 
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Fig. 6. — g vs. g - i diagrams of 100 ZAMS from Siess et al. (2000), drawn adopting increasing distance and visual extinction (left panel); and of the SDSS- 
OSIRIS sources falling in the OSIRIS field (small dots), stars with X-ray counterpart (large dots), the ZAMS from lSiess et aljj2000i) with A v = 1"' and a distance 
of 850 pc (black line), and the isochrones with age ranging from 0.5 Myrs to iOMyrs with distance and median extinction typical of cluster members (right 
panel). 




n n l I , , , , I , , , , I , , , , I , , , , I , , , , I , , , 

0.0 0.5 1.0 1,5 2.0 a. 5 3.0 



Fig. 7. — r — i vs. i - z diagram of the good - photometry stai's in the whole OSIRIS FoV (gray points). The black dots mark the colors of optical sources with 
X-ray counterparts. The black lines are the 3.5 Myr isochrone from lSiess et aljf2000l) . drawn for two different values of extinction {Ay = 2™ and Ay = 6""). 
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Fio. 8. — Histograms showing the distribution of total extinction for optical-X-ray sources in both Chandra fields and in each field alone. The vertical dotted 
lines mark the median value and the 25% and 75% quartiles. 




Fig. 9. — Left panel: gray-scale map of the median extinction for the X-ray sources wfith optical counterparts with Ay > 2"' falling in each spatial bin. Right 
panel: density map of the candidate cluster members. Stars mark the positions of the cluster members with O spectral type. 
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Fig. 10. — Color-magnitude diagrams of the good - photometry stars in the OSIRIS FOV (small gray dots). The black line is the ZAMS at a distance of 850 pc 
and an extinction Ay = 1™. The dashed lines are the isochrones for stars of 0.5, 1, 3.5, 5, and 10 Myr, at the distance but with an extinction Ay = 4™. The large 
darker dots are the optical sources with X-ray counterparts. The masses Usted on the right side of the diagram in the left panel are those predicted by a 3.5 Myr 
isochrone. 



